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Abstract 

We present a new calculation of the t cH decay rate in the standard model. We find that the 
corresponding branching fraction is B(t cH) ~ 10^^'^ 10~^^ for Mz ^ ^ 2M\y-, that is about 
6 orders of magnitude less than previously evaluated, and usually quoted in the literature. 
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Rare decays of the top quark have been extensively studied from a theoretical point of view, 
both within and beyond the standard model (SM) |l|, 0, |], ||, ^ || . The one- loop flavor-changing 
transitions, t ^ eg, t ^ cy, t ^ cZ and t — >■ cH, are particularly interesting, since they are 
sensitive to possible effects from new physics, such as supersymmetry, an extended Higgs sector 
and heavier-fermion families. In the SM, these processes are in general quite suppressed due to 
the Glashow-Iliopoulos-Maiani (GIM) mechanism [0, controlled by the light masses of the b, s, d 
quarks circulating in the loop. The corresponding branching fractions 5, = Fj/FT are further 
decreased by the large total decay width F^ of the top quark. The complete calculations of the 
one-loop flavour- changing top decays have been performed before the experimental discovery of 
the top quark in the paper by Eilam, Hewett and Soni (also based on Eilam, Haeri and Soni 
P]), where, therefore, the top mass has been taken as a parameter. Assuming rrit = 175GeV 
(according to the most recent value measured at the Tevatron rrit = (173.9 ± 5.2)GeV 0), the 
value of the total width Fj' ^ F(t ^ bW) is F^ ~ 1.55 GeV, and one gets from ref. 

B{t ^ c^) ~ 4 ■ 10"", B{t ^ C7) ~ 5 ■ 10"^^ B{t cZ) ~ 1.3 ■ 10"^^ (1) 

In the same ref. [0], a much larger branching fraction for the decay t — cH is presented as 
function of the top and Higgs masses. We show in Fig. 1 the relevant Feynman graphs for this 
channel. In particular, for nit — 175 GeV and 40 GeV ^ uih ^ 2M\y, the authors of ref. 
obtain 

S(t ^ cif) ~ 10"^ ^ 10"^ (2) 

This result was computed through the analytical formulae presented in ref. for the fourth- 
generation quark decay b' bH, in a theoretical framework assuming four flavour families. 
Such relatively large values for B{t ^ cH) look surprising, since the topology of the Feynman 
graphs for the different one-loop channels is similar, and a GIM suppression, governed by the 
down-type quark masses, is acting in all the decays. The authors of ref. argue that in the 
channel t — > cH the GIM suppression is softened by the Higgs-boson coupling to the quark. 
This statement seems puzzling, since in the t cH amplitude the relevant Yukawa coupling 
is of the order of rrit/ My/, which hardly can give rise to an enhancement of many orders of 
magnitude in the corresponding width. 

Another hint pointing to a problem with the large values of the t — >■ cH rate given above 
comes from the comparison between the rates for t ^ cZ and t cH and the corresponding 
rates for the tree- level decays t — ^ bWZ and t bWH, when mn ~ Mz- The latter can be 
considered a sort of lower-order ^^areni processes for the one-loop decays, as can be seen in Fig. 2, 
where we show the relevant Feynman graphs. In fact, the Feynman graphs for t cZ and 
t — >■ cH can be obtained by recombining the final b quark and W into a c quark in the three- 
body decays t — > bWZ and t — > bWH, respectively, and by adding analogous contributions 
where the b quark is replaced by the s and d quarks. The main point of this comparison is 
that the depletion of the t cH rate with respect to the parent t — > bWH rate is expected to 
be of the same order of magnitude of the depletion of t cZ with respect to t — > bWZ, for 
rriH — Mz- Indeed, the GIM mechanism acts in a similar way in the one-loop decays into H 
and Z. 
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An accurate determination of the t bWZ and t —>■ bWH decay rates has been carried 
out in ref. This paper also takes into account the W and Z finite-width effects that are 
crucial, since the actual value of rrit is (or can be) just near the bWZ (or bWH) threshold. For 
rriH — Mz, the two widths are comparable. In particular, for rrit ~ 175GeV, one has [jl| 

B{t bWZ) ~ 6 • 10"^ Bit bWH) ~ 3 ■ 10"^ (3) 

Assuming B{t — > cH) ~ 6 ■ 10~® that is the value taken from ref. for niH — Mz, one then 
gets for the ratio of the one-loop and tree-level decay rates 

Bit ^cH) 

th = ^ — ~ 2 (41 

- B{t ^ bWH) ^ ^ 

to be confronted with 

Bit ^ cZ) ^ 
On the other hand, it seems natural to connect the values of vh and rz to the quantity 

( ^IKM^X - 10-' (6) 



v/2' ^°'M, 
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W , 



(where Vij are the Kobayashi-Maskawa matrix elements [jTO|) arising from the higher-order in 
the weak coupling and the GIM suppression mechanism of the one-loop decay width. The large 
discrepancy between the value of the ratio in eq. (^) and what is expected from the factor 
in eq. (|^), which on the other hand is supported by the value of rz, is a further indication that 
the values for Bit cH) reported in eq. (^ can be incorrect. 

In order to clarify the situation of the t cH decay, we recomputed from scratch the 
complete anal5d;ical decay width for t — > cH (see ref. [|ll|] for details and analytical results). Our 
calculation has been done putting = 0, which automatically selects in the final amplitude 
only the dominant contribution proportional to (1 -|- 75). The corresponding numerical results 
for Bit — i> cH), when rrit = 175GeV and r(t — >■ bW) ~ 1.55 GeV, are reported in Table 1, in the 
column labeled by (m^ = 0), and in Fig. 3 (solid line). We used Mw = 80.3GeV, rrih = 5GeV, 
rUs = 0.2GeV, and for the Kobayashi-Maskawa matrix elements iV^ftVcbl = 0.04, according to 
ref. |T2[. Furthermore, we assumed |Vj*Vcs| = iV^bKfel- As a consequence, the rud dependence 



in the amplitude drops out. 

Our results are several orders of magnitude smaller than the ones reported in ref. 0. In 
particular, for ran — Mz we obtain 

BneUt ^cH) ^1.2- 10-13 (7) 

to be compared with the corresponding value presented in ref. [0] 

Boidit ^cH)c^6- 10-^ (8) 

The new value of Bit cH) in eq. (0) gives now r^ ^ 4: ■ 10^^. 
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In order to trace back the source of this inconsistency, we performed a thorough study of 
the analytical formula in eq. (3) of ref. P], for the decay width of the fourth-family down- type 
quark b' — bH, that is the basis for the numerical evaluation of B{t cH) presented in ref. [0]. 

The results of this study can be summarized in the following way. 



• We recomputed and succeeded in reproducing the complete analytical expression for 
r(5' —>■ bH) in eq. (3) of ref. |p, including the relevant form factors a and (3, govern- 
ing, respectively, the (1 -|- 75) and the (1 — 75) part of the amplitude. 

• We find good agreement with all the numerical results on the decay b' bH presented 
in ref. [|], and based on eq. (3) of the same reference. 

• We performed the numerical evaluation^] oi Bit ^ cH) by applying eq. (3) in ref. to 
the up-type quark decay t — >• cH. For rric = 0, we find complete agreement with our 
independent evaluation presented in the first column of Table 1. For completeness, in the 
last column of Table 1, we present the effect of assuming in B{t cH) a realistic value 
of the charm-quark mass rric, according to eq. (3) in p]. One can see that this effect is 
smaller than about l%o for uih ^ ISOGeV. 

• The main result of the study is that our numerical calculation disagrees by either 5 or 6 
orders of magnitude (depending on rrit and mn) with the values of B{t —>■ cH) presented 
in Fig. 3 of ref. 0. The latter are also based on eq. (3) of ref. Q. In Fig. 4jJ we present 
the updated version of Fig. 3 of ref. 

A possible explanation for this situation could be some error in the computer code used by the 
authors of ref. |^ to work out their Fig. 3. 

It is interesting to apply the complete calculation for the t cH width in the unrealistic 
parameters range 

M'^:^ml:^ml-:^ml, (9) 
where it can be compared with the results of an approximated matrix element, that is valid in 



the same parameters range |15 



^^^^ ^ Jl^^**^^^^!^^^^^^^^^ + 7')«*(P2). (10) 

This simple expression shows the basic dependence of the amplitude on both the top Yukawa 
coupling, through the factor rrit/Mw, and the remnants of the GIM cancellation, giving rise to 
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the factor V^^Vcb^-- A straightforward evaluation of the t — >• cH width according to eq. (p!0|) 
gives, for mt ^ 30GeV and mn <^ mt, 

TapproS ^ cH) ~ 4.3 ■ 10"^^GeV^, (11) 

*We checked this resuh by computing the loop scalar integrals by both our own routines and the FF library 
Fortran routines jl^ . 

^The top total decay width has been computed taking into account the complete rrib effects 
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to be compared with 

Tecactit ^ cH) ~ 6.7 ■ lO-^'^GeV. (12) 

In the real case of ~ 175GeV, there are two main hmitations of the approximation in eq. ([T0|) . 
First, the amphtude in eq. (|T^) does not take into account the absorbitive contributions arising 
when nit ~ M]y + rrit,. Secondly, it neglects the terms of order m'l/M'^ (note that both the 
terms of order mf/M^^ and the ultraviolet singularities are cancelled by the GIM mechanism). 

In Fig. 5, we compare as a function of rrit the behaviour of our calculation for r(t — »• cH) 
(solid line) with the approximated one obtained through eq. (p!OD (dashed line). Here, we 
assume the unphysical value mn = IGeV, in order to be able to explore the approximation 
down to very low rrit, maintaining the parameters in its validity region. It is worth noting 
that the simple amplitude of eq. ( p!0D reproduces very well the exact behaviour of r(t cH) 
for rrit ~ 40GeV. Note also the peak produced in the complete calculation by the opening of 
the absorbitive contribution for rrit ~ + that is obviously absent in the approximated 
estimate. In Fig. 3, we also display, for the sake of comparison, the behaviour for B(t ^ cH) 
calculated through the approximated formula eq. (p!OD and extrapolated out of its validity 
region, at the fixed value rrit = 175GeV, as a function of mn- For values of the Higgs mass 
not excluded experimentally (presently mn ~ 89GeV ||16|), the approximated amplitude gives 



results differing from the exact ones by less than a factor 4. The change in the slope of the 
exact result for rriH — 160GeV is the effect of the threshold for a new absorbitive contribution 
to the amplitude for mn ~ 2M\y. 

In conclusion, in this letter we have pointed out that the widely quoted result of ref. 
establishing a relatively large branching ratio for the decay t — > cH in the SM has been 
overestimated by several orders of magnitude. On the other hand, we agree with the analytical 
formulae in by which this result was worked out. The corrected numerical estimates are 
shown in Table 1. We find B{t ^ cH) 1 ■ 10'^^ 4 ■ 10"^^ for < ^ 2Mw, that is 
about 6 orders of magnitude less than previously evaluated, and usually quoted in the literature. 
Such a small rate will not be measurable even at the highest luminosity accelerators that are 
presently conceivable. An eventual experimental signal in the rare t decays will definitely have 
to be ascribed to some new physics effect. 

We wish to warmly thank V.A. Ilyin for discussions and suggestions. 



After submission of this paper and appearance of the preprint in the hep-ph bulletin, we 
received a message by one of the authors of ref. (J.L.H.), saying that they were aware of the 
error in B{t ^ cH), which indeed arose from a mistake in the numerical code. Unfortunately, a 
proper erratum never appeared. The correct numbers (with which we agree) have been shown 
in Fig. 21 of ref. where, however, the disagreement with the previous values in ref. has 
not been stressed. 
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Figure 1: Feynman graphs for the decay t — > cH in the unitary gauge {rric = is assumed). 




Figure 2: Feynman graphs for the decay t — > bWH {t — > bWZ). 
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mu (Gel/) 


B{t cH) {rric = 0) 


B{t^ cH) (m^ = l.bGeV) 


60 


0.2557-10-12 


0.2558 - 10-12 


70 


0.1986-10-12 


0.1986- 10-12 


80 


0.1532- 10-12 


0.1532 - 10-12 


90 


0.1169- 10-12 


0.1169 - 10-12 


100 


0.8775-10-13 


0.8777-10-13 


110 


0.6451 - 10-13 


0.6452 - 10-13 


120 


0.4605- 10-13 


0.4605 - 10-13 


130 


0.3146- 10-13 


0.3146 - 10-13 


140 


0.1999- 10-13 


0.1998 - 10-13 


150 


0.1107- 10-13 


0.1105 - 10-13 


160 


0.4424 - 10-1^ 


0.4410 - 10-1^ 



Table 1: Branching ratio for the decay t — > cH versus ttlh for both massless and massive 
charm quark. In the column (m^ = 0), we show the values obtained by our own calculation, 
in the column {m^. = 1.5GeV) we give, for comparison, the values obtained by our numerical 
evaluation of the formulae in ref. We assume rrit = 175GeV. 
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Figure 3: Branching ratio for the decay t cH versus tuh, for = 175GeV and r(t 
bW) = 1.55GeV. The sohd hue is the result of our calculation. The dashed line shows the 
approximated behaviour according to eq. (0). 



9 




10 




Figure 5: Comparison between the widths of t — cH , as functions of m^, at fixed (unreahstic) 
iriH = IGeV, and rric = 0, calculated through our complete formulae (solid line) and the 
approximate expression in eq. (|TU]) (dashed line). 
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